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I. INTRODUCTION 

Heavy ion collisions have received significant atten- 
tion in recent years. Various possible probes have been 
studied in order to detect the signatures of quark gluon 
plasma (QGP). Study of direct photon and dilepton spec- 
tra emanating from hot and dense matter formed in ultra- 
relativistic heavy ion collisions is a field of considerable 
current interest. Electromagnetic probes have been pro- 
posed to be one of the most promising tools to charac- 
terize the initial state of the collisions [ij . Because of the 
very nature of their interactions with the constituents of 
the system they tend to leave the system almost unscat- 
tered. In fact, photons (dilepton as well) can be used 
to determine the initial temperature, or equivalently the 
equilibration time. These are related to the final multi- 
plicity of the produced hadrons in relativistic heavy ion 
collisions (HIC). By comparing the initial temperature 
with the transition temperature from lattice QCD, one 
can infer whether QGP is formed or not. 

Photons are produced at various stages of the evolu- 
tion process. The initial hard scatterings (Compton and 
annihilation) of partons lead to photon production which 
we call hard photons. If quark gluon plasma (QGP) is 
produced initially, there are QGP-photons from thermal 
Compton plus annihilation processes. Photons are also 
produced from different hadronic reactions from hadronic 
matter either formed initially (no QGP scenario) or re- 
alized as a result of a phase transition (assumed to be 
first order in the present work) from QGP. In addition 
to that there is a large background of photons coming 
from tt'^ and 77" decays. If this decay contribution is sub- 
tracted from the total photon yield what is left is the 
direct (excess) photons. 

These apart, there exits another class of photon emis- 
sion process via the jet conversion mechanism (jet-plasma 
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interaction) which occurs when a high energy jet in- 
teracts with the medium constituents via annihilation 
and Compton processes. It might be noted that this 
phenomenon (for Compton process) has been illustrated 
quite some time ago [3| in the context of estimating pho- 
tons from equilibrating plasma. There, it is assumed 
that because of the larger cross-section, gluons equili- 
brate faster providing a heat bath to the incoming quark- 
jet. A comparison of the non-equilibrium photons (equiv- 
alent to photons from jet-plasma interaction) with the 
direct photons (thermal) shows that this contribution re- 
mains dominant for photons with px upto 6 GeV. How- 
ever, while evaluating jet-photon the authors in Ref. Q 
assumes that the largest contribution to photons corre- 
sponds to ^ PqiPq)- This implies that the annihilating 
quark (anti-quark) directly converts into a photon. In 
the present work, we calculate photons from jet-plasma 
interaction relaxing the above assumption and including 
the energy loss of the participating jet where the strong 
coupling constant (as) is taken as 'running' as given in 
Ref. i. 

The phenomena of jet quenching vis-a-vis energy loss 
have been studied by several authors Q taking into ac- 
count both coUisional and radiative losses. However, in 
the present work, our main concern is to examine the 
role of coUisional loss (where as is running) in the con- 
text of photon production from jet-plasma interaction 
for the following reason. The suppression of single elec- 
tron data is more than expected which led to the re- 
thinking of the importance of coUisional energy loss in the 
context of RHIC data. A substantial amount of work has 
been done to look into this issue in recent times [1, 0- 
[Tl| . Few comments about the recent developments of 
coUisional energy loss are in order here. It is argued in 
Ref. [l^ that the coUisional energy loss is approximately 
of the same order as the radiative loss. It is also shown 
by Braun et. al [3l that the coUisional energy loss in- 
creases substantially if the strong coupling is treated as 
function of temperature and momentum and if, in addi- 
tion to <-channel process, the inverse Compton reaction is 
considered. In a most rcccnt calculation, using a reduced 
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screening mass and running coupling the collisional en- 
ergy loss is six times larger than that with the constant 
coupling [l^l ■ It explains single electron Raa quite well. 
However, it fails to account for the elliptic flow, V2 of the 
electron. Effective resonance with LO-pQCD model 
also improves the collisional energy loss and the single 
electron data is well reproduced. It is also important to 
note that only the radiative energy loss fails to account 
for the single electron data at RHIC [l^. On the other 
hand, the authors of Ref. [l^ claims that the collisional 
energy loss is sub-leading. However, in order to see the 
effects of energy loss on jet-photon one should also in- 
corporate the radiative energy loss for completeness and 
this has to be done in the same formalism in a realistic 
scenario. 

Thus, it is apparent that the issue of the relative im- 
portance of the mechanism of energy loss in the context 
of RHIC data is not settled yet. We shall re- visit the 
importance of collisional energy loss in the context of 
photons from jet-plasma interactions. Moreover, in high 
temperature (T) effective field theory the coupling con- 
stant, Og, is taken to be a function of temperature only, 
which may be justified when T » Aqcd- However, in 
the case of relativistic heavy ion collisions, temperature 
is not the only scale, there is the momentum scale (k) 
also. One has to take into account the case when k ^ T 
and treat as to be function of both k and T By 
incorporating this fact it is shown that the energy loss 
is by a factor of 2 — 4 more than the case when as is 
constant. It is for this purpose we concentrate on the 
collisional energy loss 0,113 with the formalism given in 
Ref. Q to calculate photons from jet-plasma interaction. 
For completeness, we also include the radiative loss in an 
effective way. 

The organization of the paper is as follows. We give 
a brief description of jet-photon production in QGP in 
section HA. The evolution of jet quark and photon pT 
distributions are discussed in sections IIB and IIC re- 
spectively. Section HI is devoted to the discussions of 
results and finally, we summarize in section IV. 



II. FORMALISM 

A. Jet-Photon Rate 

The lowest order processes for photon emission from 
QGP are the Compton scattering {q{q) g — ?> q{q) 7) and 
annihilation {qq — 57) process. The total cross-section 
diverges in the limit i or w — >■ 0. These singularities 
have to be shielded by thermal effects in order to ob- 
tain infrared safe calculations. It has been argued in 
Ref. 



18j that the intermediate quark acquires a ther- 



mal mass in the medium, whereas the hard thermal loop 
(HTL) approach of Ref. [l^ shows that very soft modes 
are suppressed in a medium providing a natural cut- 
off kc ^ gT. We assume that the singularities can be 
shielded by the introduction of thermal masses for the 



participating partons. Apart from the thermal interac- 
tions of the plasma partons, interaction of a leading jet 
parton with the plasma was found to be a very important 
source of photons. 

The differential photon production rate for this process 
is given by: 
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where, \A4\'^ represents the spin averaged matrix element 
squared for one of those processes which contributes in 
the photon rate and M is the degeneracy factor of the 
corresponding process, fj^t , /2 and /a are the initial state 
and final state partons. /2 and /s are the Bose-Einstein 
or Fermi-Dirac distribution functions. 
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B. 



Fokker - Planck Equation: Parton transverse 
momentum spectra 



In the photon production rate (from jet-plasma inter- 
action) one of the collision partners is assumed to be 
in equilibrium and the other (the jet) is executing ran- 
dom motion in the heat bath provided by quarks (anti- 
quarks) and gluons. Furthermore, the interaction of the 
jet is dominated by small angle scattering. In such sce- 
nario the evolution of the jet phase space distribution is 
governed by Fokker-Planck (FP) equation where the col- 
lision integral is approximated by appropriately defined 
drag and diffusion coefficients. 

As mentioned already in the introduction that the 
quark jet here is not in equilibrium. Therefore the cor- 
responding distribution function (/jet) that appears in 
Eq. ([T|) is calculated by solving the FP equation. The en- 
ergy loss is represented by the drag coefficient (see later) . 
The FP equation, can be derived from Boltzmann equa- 
tion if one of the partners of the binary collisions is in 
thermal equilibrium and the collisions are dominated by 
the small angle scattering involving soft momentum ex- 
change H, [ij, l20l - [2^ . For a longitudinally expanding 
plasma, FP equation reads 0, [2^ : 
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Here v is the appropriate degeneracy factor. Note that 
the coefficient Ai is related the drag coefficient rjhy Ai ~ 
•qpi, where rj = {1/ E)dE / dx. 

Now, Bij can be decomposed into longitudinal and 
transverse components: 
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Explicit calculation shows that the off diagonal compo- 



nents of Bii vanish and we have, 
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where, Bi = d{{Ap^)^) /dt, Bt = d{{ApTY)/dt, repre- 
sent diffusion constants along parallel and perpendicu- 
lar directions of the propagating partons. Evidently, Ai 
(Btj) is infrared singular. Such divergences do not arise 
if close and distant collisions are treated separately. For 
very low momentum transfer the concept of individual 
collision breaks down and one has to take collective exci- 
tations of the plasma into account. Hence there should be 
a lower momentum cut off above which bare interactions 
might be considered. While for soft collisions medium 
modified hard thermal loop corrected propagator should 
be used 0, [M] ■ It is evident that Eq. ^ actually gives 
dE/dt or the energy loss rate that can be related to 
the drag coefficient. 

However, in the above treatment, the infra-red cut-off 
is fixed by plasma effects, where only the medium part is 
considered, completely neglecting the vacuum contribu- 
tion leading to ambiguity in the energy loss calculation. 
If the latter part is taken into account the strong coupling 
should be running. Thus for any consistent calculation 
one has to take into con sideration this fact. In that ease 
as = as{k,T) (fc = \/\uj'^ — q^\ in this case), and the 
above integrals must be evaluated numerically where the 
infra-red cut-off is fixed by Debye mass to be solved self- 
consistently: 
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We reiterate that the matrix elements in Eqs. ^ and ([7]) 
contains the strong coupling which we take as running, 
i. e. as = as(\/|w^ — <i^\i T). We chose the follow- 
ing parameterization of cts which respects the perturba- 
tive ultra-violet (UV) behavior and the 3D infra-red (IR) 
point 0: 
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FIG. 1: (Color online) pr distribution of parton at (a) RHIC 
energy {Ti = 0.446 GeV and n = 0.147 fm/c) and at (b) LHC 
energy [Ti = 0.897 GeV, n = 0.073 fm/c). 
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with k — \/\uP' — q^\ in this case. The parameters a, h 
and As are given by a = 9.07, h = 5.90 and As = 0.263 
GeV. For the limiting behavior (k << T) of the coupling 
we choose, 



ui = a*3d(l + exp{~U3)) 



(10) 



Here a*3d and a* s denote the values of the IR fixed 
point of SU (3) Yang-Mills theory in d = 3 and d = 4 
dimensions, respectively. The remaining four parame- 
ters (u2 = 5.47, M3 = 6.01, W2 = 10.13 and 1)3 = 9.27) fit 
the numerical results for pure Yang-Mills theory obtained 
from the RG equations in Ref. [27| . 

So far we have discussed about the coUisional energy 
loss which, although, dominates at lower energies, is not 
the only mechanism of energy loss. As the energy in- 
creases radiative energy loss starts to dominate and hence 
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cannot be neglected. In case of jet-photon production, 
since the photon energy is ahnost equal to the jet en- 
ergy, one has to include the radiative loss to account for 
the high px photons. However, in order to see the ef- 
fects of both the coUisional and radiative energy losses, 
one must develop a formalism in which both the mech- 
anisms can be taken into account in a consistent man- 
ner. The two mechanisms are not entirely independent, 
i.e., the collisional loss may influence the radiative loss. 
Thus both should be included to calculate transport coef- 
ficients. Since there is no rigorous way to implement this, 
the approximate way is to define effective drag (diffusion) 
in the following manner: 
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where ^ '^^ 



dx ) coll ''^^^ ^® calculated from Eq. (|4]) keeping in 
mind that rjpi = Ai. The collisional (differential) energy 
loss {dE/dx)co\\ can be calculated from Eq. Q using the 
method described in Q. For running as it is given by. 
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Where, s ~ 2ET, E is the energy of the incident quark. 
The radiative energy loss is given by. 
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where N{E) is an energy dependent factor and N{E ^ 
oo) = 4 if kinematic bounds are neglected [2^. It is 
important to point out here that N{E) = 7.3, 10.1,24.4 
for E 500, 50^5 GeV respectively and N{E ^ oo) = 4 
(see the Ref. [2g). L is the distance traversed by the 
jets in the plasma. Similarly we can define the effective 
diffusion coefHcients. 

Having known the drag and diffusion, we solve the 
FP equation using Green's function techniques: If 
P{p^t\po,ti) is a solution to Eq. ([3]) with the initial con- 
dition 
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the full solution with an arbitrary initial condition can 
be obtained as [l^] 

.f{t,p) = J d''poP{p,t\ph,U)fo{po) (15) 

where for the initial condition f{t ~ 0,p) = fo{pa) and 
P{p,t\pQ,ti) is the Green's function of the partial differ- 
ential Eq. (2). 

We assume here that the plasma expands only lon- 
gitudinally (Bjorken expansion scenario [1^). The rea- 
son is the following. The transverse expansion will have 



two effects on the parton energy loss: (i) The expanding 
geometry will increase the duration of propagation, (ii) 
the same expansion will cause the parton density to fall 
along its path. These two effects partially compensate 
each other and the energy loss is almost the same as in 
the case without transverse expansion [s^ . Since we are 
considering the central rapidity region {\r]\ < 0.35) the 
arguments given above is not applicable in the case of 
longitudinal expansion scenario. 

The solution with an arbitrary initial momentum dis- 
tribution can now be written as [l4| , 
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We use the initial parton distributions (at the forma- 
tion time ti) taken from 0, [3l|: 



dN 



KNo 



d'^Pordyt, 



\yo=o 



(1 



POT \n ' 
B I 



(17) 



where K \s a, phenomenological factor ('^1.5 — 2) which 
takes into account the higher order effects. The values of 
the parameters are listed in Table. H] We note that the 





RHIC 


LHC 
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q 
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iVo [l/GeV] 


5.0 X 10^ 


1.3 X 10^ 


1.4 X 10" 


1.4 X 10' 


P [GeV] 


1.6 


1.9 


0.61 


0.32 


a 


7.9 


8.9 


5.3 


5.2 



TABLE I: Parameters for initial parton pr distribution. 

parametric form of Eq. (jl7p may not represent the true 
picture of the jet pr distribution. In recent years, more 
sophisticated calculations have been done using different 
parameterizations of the parton distribution functions. 
The pt distribution used here might differ substantially 
from these calculations and the results presented here is 
correct upto a factor that might come from using the 
more state of the art calculation. 



C. Space time evolution 

In order to obtain the space-time integrated rate we 
first note that the phase space distribution function for 
the incor ning jet in the mid rapidity region is given by 
(see Ref. jSy for details) 
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where 



dN 



d^PTdy iPT^ can be obtained from Eq. ([TB)) . ti is 
the jet formation time and Vq is the spin-color degeneracy 



factor, zq is the jet formation position in the direction of 
QGP expansion and is the initial jet production 

probabiUty distribution at the initial radial position Wr 
in the plane zq = 0, where 
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= y^(rcos0 - t')^ + ^2 sin^(?!) for - (19) 

and (j) is the angle in the plane zq = between the di- 
rection of the photon and the position where this photon 
has been produced. We assume the plasma expands only 
longitudinally. Thus using d'^x ~ rdrdt' d(pdz and the ex- 
pression for fjet from Eq. [17] we obtain the transverse 
momentum distribution of photon as follows [32l, fssj : 
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fjet is the distribution function of the jet quark (see 



Eq. ([181)') rest of the distribution functions i.e f2, fs 
are Fcrmi-Dirac or Bose-Einstein distributions. depen- 
dence occurs only in V^Wr)- So the integration can be 
done analytically as in Ref. [s^ . The temperature profile 
is taken from Ref. [33 |. 

Besides the thermal photons from QGP and hadronic 
matter we also calculate photons from initial hard scat- 
tering from the reaction of the type hAhs 7 ^ using 
perturbative QCD. We include the transverse momentum 
broadening in the initial state partons [34il35| . The cross- 
section for this process can then be written in terms of 
elementary parton-parton cross-section multiplied by the 
partonic flux which depends on the parton distribution 
functions (PDF) for which we take CTEQ parameteriza- 
tion [3^. A phenomenological factor K is used to take 
into account the higher order effects. We also include 
photons from fragmentation process. 
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FIG. 2: (Color online) pr distribution of photons at RHIC 
energy with Ti = 0.446 GeV and n = 0.147 fm/c. The violet 
(magenta) curve denotes the photon yield from jet-plasma in- 
teraction with coUisional (collisional + radiative) energy loss. 
The blue curve corresponds to the case without any energy 
loss and the green curve represents the thermal contribution. 



the hadronic phase) where the quenching factor is the 
largest. At LHC energies (see Fig.[T}D) this factor is more 
as the temperature in this case is large compared to RHIC 
energies. It is seen from Eq. (j20p the photon pT distri- 
bution is directly proportional to the quark pT spectra. 
Thus the photon yield will be affected as we shall see in 
the followinff. 
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III. RESULTS 

To obtain the quark momentum distribution we use 
Eqs. ([3]), ([4|), ([7|), ^ and ([16|). The transverse mo- 
mentum distributions of quarks are shown in Fig. [T] for 
different times (proper) at RHIC and LHC energies re- 
spectively where the initial distributions are taken from 
Eq. ([TT]) . It is observed that the spectra are more re- 
duced as the time increases. It is generally assumed that 
quarks fragment into hadrons around Tc (the begining of 



FIG. 3: (Color online) pr distribution of photons at RHIC en- 
ergy with Ti = 0.446 GeV and n = 0.147 fm/c. The magenta 
(blue) curve denotes the photon yield from jet-plasma inter- 
action running (constant) Us- The black (green) curve corre- 
sponds to hard (thermal) photons. The orange represents the 
total photon yield compared with the Phenix measurements 
of photon data [37j . 

In order to obtain the photon pT distribution we nu- 
merically integrate Eq. (|20p using Eq. ([THj). The results 
for jet-photons for RHIC energies are plotted in Fig. [2] 
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where we have taken Ti = 446 MeV and ti = 0.147 fm/c. 
As indicated earher, the radiative energy loss starts dom- 
inating at higher energies of the jet, we include this in 
the calculation of photon px distribution. We find that 
the yield is decreased with the inclusion of both the en- 
ergy loss mechanisms as compared to the case when only 
collisional energy loss is considered. It is to be noted that 
when one considers collisional energy loss alone the yield 
with constant ag is more compared to the situation when 
running is taken into account. This is due to the fact 
that the energy loss in the later case is more [1] . On the 
other hand, with the inclusion of the radiative loss the 
yield decreases further. 
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Jet-th {with coll e-loss + running a j 
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FIG. 4: (Color online) pr distribution of jet-th photons at 
RHIC energy with Ti = 0.446 GeV and n = 0.147 fm/c. 
The red (blue) curve denotes the photon yield coming from 
jet-plasma interaction with collisional energy loss -|- running 
(constant) a^. 

In order to compare our results with high px photon 
data measured by the PHENIX collaboration [33], we 
have to evaluate the contributions to the photons from 
other sources, that might contribute in this pr range. In 
Fig. [3] the results for jet-photons corresponding to the 
RHIC energies are shown, where we have taken Ti = 446 
MeV and ti = 0.147 fm/c. The individual contribu- 
tions from hard and bremsstrahlung processes [ssj are 
also shown for comparison. These are estimated using 
the formalism given in Ref. [l^. The total yield com- 
prises of photons from jet-plasma interaction (with en- 
ergy loss), hard and bremsstrahlung processes, thermal 
Compton and annihilation processes. We also show in a 
separate plot (in Fig. ^ photons from jct-plasma inter- 
action corresponding to the cases with constant as and 
running as with collisional energy loss alone. It is ob- 
served that the spectra in the case of collisional energy 
loss with running coupling is depleted by a factor 1.7 — 2 
compared to the case where the strong coupling is con- 
stant. This is expected as the energy loss is more in the 



former case. The yield further reduces when both the 
mechanisms of energy loss are included. The total pho- 
ton yield consisting of jet-photon, photons from initial 
hard collisions, jet-fragmentation and thermal photons 
is compared with the PHENIX photon data [37|. It is 
seen that the data is well reproduced in our model (see 
Fig. [3]). To cover the uncertainties in the initial condi- 
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FIG. 5: (Color online) pt distribution of photons at RHIC 
energies with Ti = 0.350 GeV and Ti = 0.25 fm/c. The orange 
line corresponds to total photon yield from all the sources 
as in Fig. [3] with both energy losses (collisional -|- radiative) 
included in the jet-photon contribution. 

tions for a given beam energy, we consider another set of 
initial conditions at a lower temperature T, = 0.350 GeV 
and somewhat later initial time of = 0.25 fm/c. The 
yield for this set is shown in Fig. [5] We see that the data 
is reproduced reasonably well. 

In Fig. [6] we plot the px distribution of photons for 
the RHIC energy {T, = 0.236 GeV and n = 0.5 fm/c) 
for a lower value of dN/dy = 600. It is clearly visible 
from Fig. [6] that we can not explain Phenix photon data 
satisfactorily in the pr range 4 — 8 GeV. For the higher 
Pt range Phenix photon data is well reproduced. 

We also consider the high pT photon production at 
LHC energies. The contributions from various sources 
are shown in Figs. [7] where the jet-plasma contribution 
is calculated with running coupling constant (as) (con- 
sidering both the collisional and radiative energy losses). 
Since the initial temperature in this case is higher, the 
plasma lives for longer time. Thus the energy loss suf- 
fered by the parton is more. As a result, the difference 
between the cases with and without energy loss is slightly 
more than what is obtained at RHIC. It is observed that 
due to the inclusion of radiative energy loss along with 
collisional energy loss the jet photon yield is suppressed 
significantly. Also due to the inclusion of running cou- 
pling constant the jet photon yield is suppressed by a 
factor of 3 — 4. 
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FIG. 6: (Color online) pt distribution of photons at RHIC 
energy with Ti = 0.236 GeV and n = 0.5 fm/c. The magenta 
(blue) curve denotes the photon yield from jet-plasma inter- 
action running (constant) a^. The black (green) curve cor- 
responds to hard (thermal) photons. The orange represents 
the total yield compared with the Phenix measurements of 
photon data [stI ]. 
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FIG. 7: (Color online) Same as Fig.® at LHC energy {Ti 
0.897 GeV and n = 0.073 fm/c). 



IV. SUMMARY 



"We have calculated the transverse momentum distribu- 
tion of photons from jet plasma interaction with running 
coupling, i. e. vifith = as{k,T) where we have in- 
cluded both collisional and radiative energy losses. It is 
found that the assumption made in Ref. [2| while calcu- 
lating photons from jet-plasma interactions may not be 
good at LHC energies as we observe a difference is by a 
factor of 2 — 3. Using running coupling we find that the 
depletion in the photon px spectra is by a factor of 2 — 2.5 
more as compared to the case with constant coupling for 
RHIC energies. This is due to the fact that the energy 
loss (and hence drag and diffusion coefficients) is more 
by similar factor in the case of running coupling. Phenix 
photon data have been contrasted with the present cal- 
culation and the data seem to have been reproduced well 
in the low pT domain. The energy of the jet quark to 
produce photons in this range {4 < px < 14) is such 
that collisional energy loss plays important role here. It 
is shown that inclusion of radiative energy loss also de- 
scribes the data reasonable well. 

To check the sensitivity with the initial conditions, we 
consider two sets of initial conditions. In both the cases 
the data can be described quite well. This is due to the 
fact that both the initial conditions corresponds to the 
same dN/dy = 1150. 

As we validate our model through the description of 
Phenix photon data we also predict the high pT photon 
yield that might be expected in the future experiment 
at LHC. We notice that the inclusion of the radiative 
energy loss further reduces the yield at high pt- It is ob- 
served that the contribution from jet-plasma interaction 
is slightly more reduced as compared to the RHIC case 
as the initial temperature is higher at LHC. 

We do not consider transverse expansion as the en- 
ergy loss of the partons remains just about the same as 
the case without transverse expansion. Finally, we con- 
clude by noting that the role of running coupling constant 
should be explored in the context of other observables 
such as thermal photons, dileptons and so on. 
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